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The paper presents analysis of the single top quark production in PbPb collisions at the Large 
Hadron Collider at center-of-mass energy 5.5 TeV per nucleon pair. The analysis is performed 
with CompHEP and PYQUEN event generators. The neutron and proton content in the nuclei 
is taken into account. NLO precision has been implemented to simulate kinematic properties and 
rate of single top production. The modification of different characteristics of single top quark 
decay products due to interactions of jet partons in quark-gluon medium, and the specific charge 
asymmetry of top/anti-top quark yields due to the isospin effect are evaluated. 


I. INTRODUCTION 

One of the important tools to study the properties of 
hot and dense matter created in ultrarelativistic heavy 
ion collisions is a QCD jet production. Medium-induced 
energy loss of energetic partons, the so-called jet quench¬ 
ing, is supposed to be very different in cold nuclear mat¬ 
ter and in quark-gluon medium, and leads to a num¬ 
ber of phenomena which have been already seen in the 
RHIC data on high momentum particle production in 
gold-gold collisions (see, e.g., reviews m and references 
therein). The LHC heavy ion program makes it possible 
to probe the new frontiers of the high temperature Quan¬ 
tum Chromodynamics (QCD) with the increasing role of 
hard and massive particle production processes. The ef¬ 
fective reconstruction of energetic jets and other hard 
probes (like W-boson and b-tagged jets) at LHC ener¬ 
gies becomes possible even in heavy ion environment. A 
number of interesting LHC results from PbPb runs at 
V^nn = 2.76 TeV, have been published by ALICE, AT¬ 
LAS and CMS collaborations (the overview of the results 
from the first year of heavy ion physics at LHC can be 
found in [6]). In particular, the first direct observation 
of jet quenching has been done measuring the transverse 
energy imbalance in dijet CHS] and photon+jet m chan¬ 
nels. Then jet quenching effect has been manifested by a 
number of specific phenomena, such as the suppression of 
inclusive mm and semi-inclusive [T4] jet rates in cen¬ 
tral PbPb collisions compared to peripheral events and 
proton-proton interactions, similar suppression for jets 
from b-quark fragmentation m , azimuthal anisotropy of 
jet spectra m , and medium-modifed jet internal struc¬ 
ture mm (see P3 for some overview). At that a num¬ 
ber of theoretical calculations and Monte-Carlo simula¬ 
tions in many different approaches were attempted to re¬ 
produce the jet quenching observables in PbPb collisions 
at the LHC, and to extract by such a way new infor¬ 
mation about created medium and partonic energy loss 
mechanisms j2QN4l] . The review of the current status of 
jet quenching theory can be found e.g. in [42]. 

Single top quark production is a new interesting hard 
probe of quark gluon matter at the LHC. Top quarks are 
produced and decay at very early stage of the nuclear re¬ 
action before the dense matter formation. Then partonic 
products of top quark decay loose the energy propagating 


through the QCD-medium, and so their different kine¬ 
matic distributions get modified. In particular, smearing 
and decreasing mean and maximum values of invariant 
mass distributions of three jets from top quark decay 
(one b-jet and two jets from W-boson) and dijets from 
W-boson decay in PbPb collisions as compared with pp 
interactions are predicted for top anti-top pair produc¬ 
tion [43]. Single and pair top-quark production cross 
sections in proton-lead and lead-lead collisions at the en¬ 
ergies of LHC and Future Circular Collider (FCC) have 
been estimated recently in [44] with next-to-leading-order 
perturbative QCD calculations including nuclear parton 
distribution functions. 

In our paper we analyze the influence of nuclear ef¬ 
fects on single top decay pattern in PbPb collisions at 
nominal LHC energy V^nn = 5.5 TeV in the frameworks 
of CompHEP [47] and PYQUEN [48] event generators. 
The cross section calculation of single top production 
process is done with MCFM [49], subsequent generation 
of nucleon-nucleon events is done with CompHEP, while 
PYQUEN being used to simulate in-medium rescattering 
and energy loss of top quark partonic decay products. 
The neutron and proton content in Pb nuclei, and nu¬ 
clear parton distribution functions have been taken into 
account. The medium-induced modification of different 
characteristics of single top quark decay products (in¬ 
variant mass distribution of W-boson and b-jet from top 
quark decays, the fraction of top quarks observed with 
associated partner jet, transverse momentum distribu¬ 
tion of associated jets), and the charge asymmetry of 
top/anti-top quark yields due to the isospin effect are 
investigated. 


II. SIMULATION OF SINGLE TOP QUARK 
PRODUCTION 

At hadron and lepton colliders top quarks are produced 
either inpairs or singly. The representative diagrams for 
the single top production are shown in Fig. []] 

Three mechanisms of the single top production are dis¬ 
tinguished by the virtuality of the W-boson involved: 
t-channel (Q^ < 0), 5-channel (Q^ > 0), associated 
tW (Qw = Mw)- At LHC energies, most of single top 
quarks are produced with associated partner jet in t 
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FIG. 1. The representative diagrams for the single top pro¬ 
duction at hadron colliders 
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TABLE I. For all scenarios general simulation procedure 
(CompHEP and MCFM with MSTW2008NNLO + PYTHIA) 
is implemented. Simulation procedure differences are shown 
in the tabel. 


channel. In our studies we are taking into account only 
^-channel, as the most significant one. 

We consider three scenarios: (1) without any nuclear 
effect; ( 2 ) with initial state effect (nuclear shadowing); 
(3) with both initial and final (jet quenching) state ef¬ 
fects. This notations are used on all the plots and are 
described in details in this section and Table [H 


A. Initial state 

We use CompHEP for effective NLO generation [50] of 
single top events at v^nn = 5.5 TeV, where exact NLO 
cross sections are obtained with MCFM [39]. CompHEP 
is a general, tree level generator, which allows one to 
study various SM and BSM processes 2 —>> N (up to 
N = 6 ) in the framework of the usual technique of Feyn¬ 
man diagrams squared for different models. It gener¬ 
ates, squares and symbolically calculates a set of Feyn¬ 
man diagrams for a given process and creates a numeric 
Monte-Carlo generator for the process. This generator 
allows one to compute cross sections (with applied cuts), 
to build distributions and to generate events with partons 
in the final state; the initial partons are convoluted with 
parton distribution functions (PDF). In nucleus-nucleus 
A A collisions, the parton flux is enhanced by the num¬ 
ber A 2 of nucleons in the nucleus and modulo small 
(anti)shadowing effects in the nuclear PDF. In our cal¬ 
culations we have used MSTW2008 PDF |5T] at NNLO 
with nuclear corrections PDF EPS09 .52] from LHAPDF 



FIG. 2. Experimental data for charge asymmetry 
(Nw+ — Nw-)/(Nw+ + Nw-) is shown for pp(blue open 
squares) and PbPb(red-filled circles) at yG = 2.76 [55]. Mod¬ 
eling for pp collisions, first scenario, is shown with blue solid 
line. Second, PbPb collisions - green dashed line. The third 
one, PbPb collisions with final state effects are shown with 
red dotted line. Scenario details are described in Tab. [T| 


library [53]- It allowed to take into account nuclear 
(anti)shadowing modifications of the bound relative to 
free nucleons as well as the different isospin (u- and d- 
quark) content of Pb ion given by its different proton 
(Z=82) and neutron (N=126) numbers. The latter effect 
is important in the case of electroweak isospin-sensitive 
processes like W, Z, prompt photon and single top quark 
production. 

The parton showering and hadronization for generated 
events have been performed with PYTHIA.6.4 event gen¬ 
erator [54] . The correctness of our simulation procedure 
Tab. [T] is confirmed in a good agreement with the CMS 
data on W-production cross section and rapidity depen¬ 
dence of W + ,W~ charge asymmetry (Fig. [ 2 | in PbPb 
collisions at y/s NN = 2-76 TeV [55] . 

We consider the channel with leptonic (t Wb -A bvi , 
l = fi. e) top quark decay modes. The possibility to re¬ 
construct W-bosons in PbPb collisions with CMS detec¬ 
tor using missing- Et technique has been revealed in m- 
The efficiency of b-tagging in PbPb collisions is estimated 
on the level of rsj 50% US]. Thus the reconstruction of 
top quark in PbPb collisions in the leptonic decay mode 
looks quite realistic. The reconstruction of hadronically 
decayed single top quarks (t Wb bqq') seems prac¬ 
tically impossible even in pp collisions due to huge QCD 
background. 

In order to estimate the expected event rate for re¬ 
alistic geometrical acceptance and kinematic cuts, the 
selection requires pseudorapidity coverage | 77 |< 5 for 
light-flavor quark jets, presence of at least one b-tagged 
jet with | 77 | < 2.5 and | 77 |< 2.5 for leptons (muons 
and electrons), where 77 = — ln[tan(0/2)]. The criteria 
corresponds to the conditions of ATLAS and CMS ex- 
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FIG. 3. Distribution of invariant top mass reconstructed from 
b-jet and W-boson for different scenarios. Scenario details are 
described in Tab. [T| 

periments. The cuts on transverse momenta of leptons 
p l T > 20 GeV/c and jet transverse energy > 20 GeV 
were applied. 

For simplicity jet clustering is done with PYCELL 
Pythia routine. Cone based criteria A R < 0.4, where 
A R = y/ A rj 2 + A(p 2 , is used to associate jet to parton. 
In single top quark production topology presupposes two 
hard jets with high Pj-, b-jet from top quark decay and 
light quark jet and one soft jet with low Pt from addi¬ 
tional b-quark. 

We estimate efficiency of kinematic acceptance as 
^^7 ’ where N pass is number of events passed kinematic 
cuts, Ntotal is full number of events generated by Monte- 
Carlo, and it is ^ 65%. B-tagging acceptance is esti¬ 
mated as 

e X Nib jet + (1 — (1 — e ) ) x N2b jet 
Nib jet + N2bjet 

where Nibjet, N 2 bjet are number of events passed kine¬ 
matic cuts with reconstructed one and two b-jets cor¬ 
respondingly, e is b-tagging efficiency. For e = 0.5 we 
can estimate acceptance ^ 52%. Total acceptance as 
- 31 - 35%. 

Single top t-channel production NLO cross section is 
calculated with MCFM and is 38.6 pb. We assume 
equal branching | for muon and electron and total ac¬ 
ceptance as 35%. Then the estimated pp cross section 
for single top quark production in leptonic mode (elec¬ 
tron and muon channels without any contribution from 
tau-channel) satisfying above cuts is ~ 3 pb, and PbPb 
cross section is estimated as ~ 0.13/ib. Our numbers 
are compatible with ones obtained in [44] (taking into 
account that different nucleon PDF and little different 
kinematic cuts were used). The corresponding event 
rate in a one month PbPb run (T = 10 6 seconds), with 
nominal luminosity L = 10 27 cm -2 s -1 , is estimated as 
N ev = TcrpbPb^ ~ 130 (this corresponds to lnb -1 inte¬ 


grated luminosity data). In the case of increasing LHC 
luminosity for lead beams in the future (upgrade plans 
for HL-LHC) by a factor of 10, the number of detected 
events with single top quark production may increase up 
to ~ 1300. 

The main background sources for single top produc¬ 
tion in pp collisions are W+jets, tt pairs and QCD 
jets [[45], [46]]. It is natural to expect that these pro¬ 
cesses will still form the main background for single top 
production in heavy ion collisions. However their sig¬ 
nificance can be modified, as well as influence of initial 
and final state nuclear effects on them may differ. Back¬ 
ground modelling and estimation of signal statistical sig¬ 
nificance is complex and important task, which however 
goes beyond the current paper topic. It could become 
the subject of our future studies. 


B. Medium-induced final state interactions 

After specifying initial partonic state from top quark 
decays with CompHEP and PYTHIA, event-by-event 
Monte-Carlo simulation of medium-induced rescatter¬ 
ing and energy loss of jet partons is performed with 
PYQUEN event generator [48]. This model is con¬ 
structed as a modification of an event obtained with the 
PYTHIA. The approach describing the multiple scat¬ 
tering of hard partons is based on accumulated energy 
loss via gluon radiation which is associated with each 
parton scattering in expanding quark-gluon fluid. It 
also includes the interference effect in gluon emission 
with a finite formation time using the modified radia¬ 
tion spectrum as a function of the decreasing tempera¬ 
ture. The model takes into account radiative and colli- 
sional energy loss of hard partons in longitudinally ex¬ 
panding quark-gluon fluid, as well as the realistic nuclear 
geometry. The event-by-event simulation procedure in 
PYQUEN includes the generation of initial parton spec¬ 
tra with PYTHIA and production vertexes at the given 
impact parameter of AA collision; the rescattering-by¬ 
rescattering simulation of the parton path in a dense nu¬ 
clear overlapping zone inducing its radiative and colli- 
sional energy loss; and the final hadronization according 
to the Lund string model for hard partons and in-medium 
emitted gluons. 

The strength of partonic energy loss in PYQUEN is 
determined mainly by the initial maximal temperature 
T 0 max of hot fireball in central PbPb collisions, which is 
achieved in the center of nuclear overlapping area at mid¬ 
rapidity. The transverse energy density in each point 
inside the nuclear overlapping zone is supposed to be 
proportional to the impact-parameter dependent prod¬ 
uct of two nuclear thickness functions Ta in this point, 
£(t*i, 7*2) oc Ta{ti)Ta{t2) (here 7*1,2 are the transverse dis¬ 
tances between the centres of colliding nuclei and the par- 
ton production vertex). The rapidity dependent spread¬ 
ing of the initial energy density around mid-rapidity 
y = 0 is taken in the Gaussian-like form. The radia- 
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FIG. 4. Distribution of transverse momentum of b-jet from top quark decay for different scenarios. Left figure is in linear scale 
right in logarithmic. Scenario details are described in Tab. [T| 
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FIG. 5. Distribution of transverse momentum of jet from light quark for different scenarios. Left figure is in linear scale right 
in logarithmic. Scenario details are described in Tab. |T] 


tive energy loss in the model roughly are proportional 
to Tq , and collisional loss roughly are proportional to 
Tq . Such strong dependence of jet quehcning effect on 
an initial temperature allows us to fix this model parame¬ 
ter within the experimental constraints on jet quenching 
data. The partonic energy loss in the model depends 
also on the proper time To of matter formation and the 
number Nf of active flavors in the medium. The vari¬ 
ation of To value within its reasonable range has rather 
moderate influence on the strength of partonic energy 
loss. The jet quenching gets stronger at larger To due 
to slower medium cooling, which implies the jet partons 
spending more time in the hottest regions, and as a re¬ 
sult the rescattering intensity gets stronger [24] . Increas¬ 
ing Nf result in larger energy loss due to the extension 
of medium density in this case. The parameter values 
T 0 max = 1.1 GeV, t 0 = 0.1 fm/c and N f = 0 (gluon- 


dominated plasma), and also the “wide-angle” paramet- 
rization of medium-induced gluon radiation were used for 
current investigation. PYQUEN with such settings (just 
with T 0 max = 1 GeV) reproduces the LHC data on the 
dijet asymmetry [24], charged particle nuclear modifica¬ 
tion factor (up to px ^ 100 GeV/c) [56], jet fragmen¬ 
tation function and jet shapes [41] in PbPb collisions at 
y^nn = 2.76 TeV. 


III. RESULTS 

Figure [3] shows the invariant mass distributions of W- 
boson and b-jet from top quark decays in minimum bias 
(integrated over all impact parameters) PbPb interac¬ 
tions and pp collisions at v^nn = 5.5 TeV. The smearing 
and decreasing mean and maximum values of this distri- 
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FIG. 6. Distribution of pseudorapidity of jet from light quark. 
Scenario notations are similar to other figures, details are de¬ 
scribed in Tab. HJ 


FIG. 8. Top/anti-top ratio Nt/Nf. Scenario details are de¬ 
scribed in Tab. [I] 



FIG. 7. Cosine of angle between lepton from top decay and 
jet from light quark associated with the top in top quark rest 
frame. Scenario details are described in Tab. [I] 


bution from top quark mass are shown in PbPb events 
relatively to pp collisions. Transverse momentum spec¬ 
tra of jets associated with top quark is presented in Fig¬ 
ure [4] Significant softening of this distribution and peak 
shifting from the half mass of W-boson is seen for PbPb 
collisions. Softening of jet from light quark is also ob¬ 
served in Figure^] and its pseudorapidity modification is 
seen on Figure [6] 

One expects that any lepton kinematic distributions 
are not practically affected by the nuclear medium. How¬ 
ever, the distribution on the angle taken in the top rest 
frame between the lepton from the top decay and the 
momentum of the light jet produced in association with 
the top is slightly affected due to jet quenching of all 
the jets not only from the top decay. The distribution 
on this important variable for separating the single top 
signal from backgrounds is shown in the Figure [7] 

Another interesting (initial state) nuclear effect, which 
is not directly related to quark-gluon matter formation 


- (i) pp 

-- (2) PbPb 

. (3) PbPb + PYQUEN 


1 1.5 

h(lep)| 


FIG. 9. Charge asymmetry ( N t — iVf)/( N t + Nf) for three 
scenarios. Scenario details are described in Tab. |U 


but reflects the neutron and proton content in Pb nuclei, 
is significant decrease in ratio of top to anti-top quarks 
cross sections in PbPb with respect to pp collisions at the 
same beam energy per nucleon pair. The predicted ratios 
of t/t yields for minimum bias PbPb events and pp col¬ 
lisions at v^nn = 5.5 TeV are ^ 1 and ~ 2 respectively. 
Dependence of these ratios from lepton r] is shown on the 
Figure [8] The same effect leads to the disappearance of 
charge asymmetry in PbPb collisions (Fig. [9|. 


IV. CONCLUSIONS 

The t-channel single top production in PbPb collisions 
at nominal LHC energy y/s NN = 5.5 TeV has been ana¬ 
lyzed in the frameworks of MCFM(calculation of cross 
sections), CompHEP (generation of events), PYTHIA 
(parton showering and hadronization) and PYQUEN 
(medium-induced partonic rescattering and energy loss) 
Monte-Carlo models. The neutron and proton content in 
Pb nuclei, and nuclear parton distribution functions have 
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been taken into account. Single top quark production has 
a large enough cross section and visible event rate for the 
nominal LHC luminosities and is open to study in PbPb 
collisions. 

Medium-induced partonic energy loss can result in 
the modification of different characteristics of single top 
quark decay products in PbPb collisions as compared 
with pp interactions. We predict smearing and decreas¬ 
ing mean and maximum values of the invariant mass dis¬ 
tribution of W-boson and b-jet from top quark decays, 
significant softening the transverse momentum spectrum 
of jets associated with top quark and shifting its peak 
from the half mass of W-boson. 

Distribution on the angle taken in the top rest frame 
between the lepton from the top decay and the momen¬ 
tum of jet from light quark produced in association with 


the top is only slightly affected and may be used for 
separation single top signal from backgrounds in future 
analysis. 

Besides due to the nuclear isospin effect, the ratio of 
top to anti-top quark cross sections becomes close to 
unity in PbPb being by a factor ~ 2 smaller than in 
pp collisions. 
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